Oscillating Universe in Hofava-Lifshitz Gravity 
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We study the dynamics of isotropic and homogeneous universes in the generalized Hofava-Lifshitz 
gravity, and classify all possible evolutions of vacuum spacetime. In the case without the detailed 
balance condition, we find a variety of phase structures of vacuum spacetimes depending on the 
coupling constants as well as the spatial curvature K and a cosmological constant A. A bounce 
universe solution is obtained for A > 0, K — ±1 or A = 0, K = — 1, while an oscillation spacetime is 
found for A > 0, K = 1, or A < 0, K — ±1. We also propose a quantum tunneling scenario from an 
oscillating spacetime to an inflationary universe, resulting in a macroscopic cyclic universe. 



PACS numbers: 04.60.-m, 98.80.Cq, 98.80.-k 



I. INTRODUCTION 

Recently Hofava proposed a power-counting renor- 
malizable theory of gravity [l|, which has attracted 
much attention over the past year. In Hofava's theory, 
Lorentz symmetry is broken and it exhibits a Lifshitz-likc 
anisotropic scaling in the ultraviolet (UV), t —> £ z t, x —> 
£x, with the dynamical critical exponent z = 3. (For this 
reason the theory is called Hofava-Lifshitz (HL) gravity.) 
It is then natural to expect that the UV behavior of the 
theory would give rise to new scenarios of cosmology @- 
0]. Earlier works have indeed revealed some interesting 
aspects of HL cosmology such as dark matter as integra- 
tion constant [j|, the generation of chiral gravitational 
waves from inflation Q , scale- invariant fluctuations with- 
out inflation [3] , and possible dark energy scenario 
There are also some discussion closely related to obser- 
vational cosmology and astrophysics such as cosmologi- 
cal perturbation [^. [Tl| - |22j . observational constraints (23r - 
l25l | , primordial magnetic field without inflation [26[ , and 
a relativistic star (23, [H| . 

Though the viability of HL gravity is still under in- 
tense debate d, [29h4jT | , we give the theory the benefit of 
the doubt and will furthermore pursue consequences of 
Hofava's intriguing idea. 

Here we focus on the dynamics of Friedmann-Lcmaitre- 
Robertson- Walker (FLRW) universe in HL gravity, which 
may provide us new aspect of the early universe. We 
study the FLRW spacetime with non-zero spatial curva- 
ture in the context of HL gravity. The non-trivial cosmo- 
logical evolution is brought by the various terms in the 
potential which are constructed from the spatial curva- 
ture 'K.ij. In particular, we find non-singular behavior in 
high curvature region, which may lead to avoidance of a 
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big bang initial singularity [J] . Consequently, many stud- 
ies on the dynamics of the FLRW universe in HL gravity 
come out within the last few years 

In the original HL gravity, the so-called detailed bal- 
ance condition is assumed. However, this condition can 
be loosened to have arbitrary coupling constants [r3l| . 
Hence both models have been so far discussed in the anal- 
ysis of the FLRW universe. As for matter components, 
a perfect fluid with the equation of state P = wp and a 
scalar field have been discussed. As a result, we can clas- 
sify the analyzed models into four types: (1) The model 
with the detailed balance condition and with a perfect 
fluid [134511, (2) The model with the detailed balance 
condition and with a scalar field 

mm ( 3 ) The 

model 

without the detailed balance condition and with a per- 
fect fluid HI, and (4) The model without the 
detailed balance condition and with a scalar field p| [59| . 

So far many works have been done assuming the de- 
tailed balance condition, and they reveal the possibility 
of singularity avoidance such as a bounce universe [3- 
Ipl l43l IH, [53 and an oscillating spacetime [4a - 
48(| . The initial singularity is avoided because of 
"dark" radiation, which is a negative a -4 term. It comes 
from higher curvature terms. The "dark" radiation was 
first introduced in the context of a brane world [62| . Al- 
though such an effect is very interesting and important, 
the "dark" radiation term may fail to avoid a singularity 
if one include radiation or massless field. The conven- 
tional radiation behaves as a -4 with positive coefficient. 
If we have a sufficient amount of real radiation, the uni- 
verse will inevitably collapse to a big-crunch singularity. 
Furthermore, if we assume that radiation field has also 
the same scaling law as gravity in the UV limit, the en- 
ergy density of radiation field changes as a -6 [23j , which 
is the same scaling law of stiff matter in the conventional 
theory. Inclusion of the positive a -6 term will kill the 
possibility of singularity avoidance by "dark" radiation. 
In order to save the present mechanism for singularity 
avoidance, one needs a negative a~ 6 term, which may be 
obtained in the generalized HL gravity model [6lj . 
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Recently some papers have discussed the case without 
the detailed balance and studied a singularity avoidance 
(a bounce universe or an oscillating behavior): One is 
by use of a phase space analysis [53l |59| , and the other is 
the case with perfect fluid with time-evolving equation of 
state (BH fool ] . The former analysis was not properly per- 
formed because they introduce the dynamical variables 
more than the degrees of freedom. In the latter case, 
although they discuss some interesting transitions, the 
assumption of the equation of state is not so clear. 

In the present paper, since there has so far not been a 
systematic and substantial analysis in cosmology based 
on this most general potential without the detailed bal- 
ance condition, we provide a complete classification of 
the cosmological dynamics. We do not include any mat- 
ter fields not only for simplicity but also to avoid unclear 
assumption. It is just straightforward to include perfect 
fluid with the equation of state P = wp (w;=constant). In 
particular, our analysis includes matter fluid with radia- 
tion and stiff matter as it is. We clarify which conditions 
should be satisfied for singularity avoidance. We also 
propose some possible scenario for a cyclic universe, i.e., 
the oscillating spacctime will transit by quantum tunnel- 
ing to an inflationary phase, resulting in a cyclic universe 
after reheating. 

The paper is organized as follows. After giving the 
generalized model of Hofava-Lifshitz gravity in § (TT1 we 
study the isotropic and homogeneous vacuum spacetime 
in § 11111 We find a variety of phase structures including a 
bounce universe and an oscillating universe. We then in- 
voke a more realistic cosmological model which may lead 
to a macroscopic cyclic universe via quantum tunneling 
from an oscillating universe. 



II. HORAVA-LIFSHITZ GRAVITY AND THE 
COUPLING CONSTANTS 

The basic variables in HL gravity are the lapse func- 
tion, N, the shift vector, Ni, and the spatial metric, gtj. 
These variables are subject to the action [l|, [6l[ 

S = 2^ / dtd3x Vd N {&K - *H L [<7ij]) , (2.1) 
where k 2 = 1 /Mp L and the kinetic term is given by 



%K = X.i i X. i i -XX 2 (2.2) 



with 



1 

2N 



(9ij - ViNj - VjJVi) 



(2-3) 



being the extrinsic curvature. The potential term T%l 
will be defined shortly. In general relativity we have 
A = 1, only for which the kinetic term is invariant under 
general coordinate transformations. In HL gravity, how- 
ever, Lorentz symmetry is broken in exchange for renor- 
malizability and the symmetry of the theory is invariance 



under the foliation-preserving diffcomorphism transfor- 
mations, 



(2.4) 

As implied by the symmetry (|2.4p it is most natural to 
consider the projectable version of HL gravity, for which 
the lapse function is dependent only on t: N = N(t) [l|. 
Since the Hamiltonian constraint is derived from the vari- 
ation with respect to the lapse function, in the pro- 
jectable version of the theory the resultant constraint 
equation is not imposed locally at each point in space, 
but rather is an integration over the whole space. In the 
cosmological setting, the projectability condition results 
in an additional dust-like component in the Friedmann 
equation [see Eq. (|3.2p below] 

The most general form of the potential is given 
by M 



2K + gi %. 

+k 2 (g 2 3L 2 + '/r/v/,'.'v.' ) + '• '.'/ i f '**' ■K^VrK.'. 



(2.5) 



where A is a cosmological constant, and %_ are the 
Ricci and scalar curvatures of the 3- metric gij, respec- 
tively, and <?i's (i = 1,...,9) are the dimensionless cou- 
pling constants. (See Appendix A for some conditions 
on these coupling constants.) 

In the original proposal [l[ Hofava assumed the 
detailed balance condition, by which the potential 
term (|2.5I) is simplified to some extent. The potential 
under the detailed balance condition is given by 



% B = - 



2(3A-1) 2(3A-1) 
(4A-1>W kV 

8(3A 
2k?\i 



2 



2k, 2 



■CuC 



where 



e ikt V k [ * 



(2.6) 



(2.7) 



is the Cotton tensor, and A^, /i and ui are constants. The 
potential (|2.6[) is therefore reproduced by identifying 



(2-8) 
(2-9) 
(2-10) 



A = 


3(3A-1) 




2 I j 2 k 2 ' 




.9i = 


-1, 






(4A-1) 2 


2 2 2 


.92 = 


4(3A-1) M 


k , g 3 = fi k , 




4^iK 2 


2k 2 10k 2 


.94 = 


o ) 95 : 










12k 2 


Zk 2 4k 2 


.97 = 


w4 , .98 


2^4' " ~ ~ ' 



(2-11) 
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and Aiy = —(3 A — l)/(/i K ). In the detailed balance 
case /i and ui are two free parameters. 

In what follows, we adopt the unit of k 2 = l(Mp L = 1) 
for brevity. 



III. FLRW UNIVERSE IN HORAVA-LIFSHITZ 
GRAVITY 

We discuss an isotropic and homogeneous vacuum uni- 
verse in Hofava-Lifshitz gravity. Note that such a vac- 



uum spacetime is not realized in general relativity. We 
will extend our analysis to anisotropic spacetime (Bianchi 
cosmology) in the separate paper. 

Assuming a FLRW spacetime, which metric is given 

by 



ds 2 



-dt 



2 ' a 2 



dr 2 



1 - Kr 2 



r 2 dn 2 



(3.1) 



with K = or ±1. We find the Friedmann equation as 



J 



K 



(3A-l)a 2 3(3A-1) 



A 



r 



9d 

n3 



(3.2) 



where H = a/a, 



9d 
9i 
9b 



= 8C. 

= 6(g 3 + 3g 2 )K 2 , 

= 12(9g 5 + 3 56 + g 7 )K 3 . 



(3.3) 



A constant C may appear from the projectability condi- 
tion and could be "dark matter" Q . For a flat universe 
[K = 0), the higher curvature terms do not give any 
contribution, and then the dynamics is almost trivial. 
Hence, in this paper, we discuss only non-flat universe 
(K = ±1). 

If A = 1, we find a usual Friedmann equation for an 
isotropic and homogeneous universe in GR with a cos- 
mological constant, dust, radiation and stiff matter. If 
gd,<?r 5 an d g s are non- negative, such a spacetime gives 
a conventional FLRW universe model. However, since 
those coefficients come from higher curvature terms, their 
positivity is not guaranteed. Rather some of them could 
be negative. As a result, we find an unconventional cos- 
mological scenario, which we shall discuss here. In what 
follows, we assume that A > 1/3, but do not fix it to be 
unity. 

In this paper, we assume C = just for simplicity. The 
Friedmann equation is written as 



2°' 



+ u{a) = 0, 



where 



iz(a) 



1 



3A- 1 



K - - 



3a 2 



9s 
3a 4 



(3.4) 



(3.5) 



Since the scale factor a changes as a particle with zero 
energy in this "potential" u, the condition u{a) < 
gives the possible range of a when the universe evolves. 
So we can classify the "motion" of the universe by the 
signs of K and A, and by the values of g r and g s . Note 



that in the case with the detailed balance condition, we 
have 



3 M 2 



th = 6( 53 + 3. 92 ) = 2(3A _ i; 

g s = 12(9.95 + 3(7 6 + g 7 )K = 



< for A > 1/3 
(3-6) 



It is some special case of our analysis, although its dy- 
namics will be completely different from generic cases 
because g s vanishes. 

We find mainly the following four types of the FLRW 
universe: 

(1) [<B<B =>$c]: Suppose u(a) < for a e (0, ar], and 
the equality is true only when a = ar- A spacetime 
starts from a big bang (s®) and expands, but it even- 
tually turns around at a = ay to contract, finding a big 
crunch (®c). is a scale factor when the universe turns 
around from expansion to contraction. 

(2) [«« =>■ oo or oo =>■ ac]: If u{a) < for any posi- 
tive values of a, a spacetime starts from a big bang and 
expands forever, or its time reversal (A spacetime con- 
tracts to a big crunch). As for the asymptotic spacetime, 
we find a oc t [Milne -m\ (K = — 1) for A = 0, while 
a cx exp( v /A/3t) [de Sitter: ds) for A > 0. We denote 
them as =>■ rW, and as ds, respectively. For the 
contracting cases, we describe them as M <BC, and 
ds => «C, respectively. 

(3) [(Bounce]: If ll{a) < for a e [ar, oo) and the equality 
holds only when a = ax, a spacetime initially contracts 
from an infinite scale, and it eventually turns around at 
a finite scale ar, and expands forever. The asymptotic 
spacetimes are the same as the case (2): M , and ds. 

(4) [oscillation]: If u (a) < for a € [a m in, Omax] and the 
equality holds only when a = a min and a = a max , a 
spacetime oscillates between two finite scale factors. 

For some specific values (or specific relations) of g r and 
<7 S , which divides two different phases of spacetimes, we 



4 



find a static universe (s ) : 

(5) [s]i A spacetime is static with a constant scale factor 
as, if u(a s ) — and u'{as) = 0. 

There are two types of static universes: one is stable 
(s B ) and the other is unstable (s u ). When we have an un- 
stable static universe, we also find the following types of 
dynamical universes with a static spacetime as an asymp- 
totic state as well: 

(6) [s u =>■ oo or oo =>5 U ]: If u(a) < for a G [as, oo) and 
the equality holds only at as, a spacetime starts from a 
static state in the infinite past, and expands forever, or it 
initially contracts from an infinite scale, and eventually 
reach a static state in the infinite future. We then have 
S u => ds, M or is, <M =>5 U 

(7) Su, or s u => «c]: If v(a) < for a G (0,as] and 
the equality holds only at as, A spacetime starts from a 
big bang and expands to a static state with a finite scale 
as, or its time reversal (A spacetime contracts from a 
static state to a big crunch). 

(8) [s u => "Bounce S u ]: If U (o) < for a G [a s ,a T ] (or 
a G [gst, as]) and the equality holds only at as and ot, a 
spacetime starts from a static state in the infinite past, 
and expands (or contracts). It eventually bounces at a 
finite scale Or, and then reach a static state again in the 
infinite future. 

For the case of A ^ 0, introducing the curvature scale 
I which is defined by 



A 
3~ 



(3.7) 



where e = ±1, we can rescale the variables and rewrite 
the "potential" 11 by the rescaled variables as 



to a 2 and a 2 = a| is the double root, we have the third 
root, which is given by 



(it 



a 



1 r 
3e 



KT2VT 



, (3.11) 



where the universe turns around (or bounces). To exist 
such a point, it must be real and positive. 
If A = 0, we find 



3s 



K a 



(3.12) 



which is found from (|3.9j) in the limit of e = 0. The cor- 
responding curve on the g r -g s plane is denoted by Tq,k- 
The radius is given by 



as 



,[0,K] 



(3.13) 



assuming Kg I > 0. 

Note that our classification depends just on g x and g s 
(or g T and g s ), apart from K and A. Since g r and <7 S are 
given by g 2 ,g3,95,96 and g 7 , but do not include g 4 ,g 8 
and gg , the fate of the universe is classified only by the 
conditions on the coupling constants of higher curvature 
terms but not on those of their derivatives such as V 

In the case with the detailed balance conditions, we 
find A < from Eq. (|2.8[) . and then obtain from Eq. 

MM 



g T = -9/4 , g s 



0. 



(3.14) 



Now we shall discuss what kind of spacetimes are real- 
ized under which conditions in the following three cases 
separately [A. A = 0, B. A > 0, and C. A < 0]. 



u{a) 



1 



3A- 1 



K 



ea 2 



3a 2 



9s 
3a 4 



(3.8) 



A = 



where a = a/i, g r = g^jt 2 , and g s = g s /( 4 - Using this 
potential and variables, we can discuss the fate of the 
universe without specifying the value of A. 

A static universe will appear if we find a solution a = 
as(> 0) which satisfies u(as) = and 11' {as) = 0. If 
A^0(e = ±l), it happens if there is a relation between 
<7 r and g s , which is defined by 



(J* 



1 

9? 



2K - 3eKg, ± 2(1 - eg r ) 



(3.9) 



This gives the curve T €t K(±) on the g T -g s plane, which 
gives the boundary between two different phases of space- 
time. The radius of a static universe is given by 



(is 



£5r 



(3.10) 



if it is real and positive. Here ± correspond to the curves 
T*e,(±)- Since u(a) = is the cubic equation with respect 



If a cosmological constant is absent, the "potential" is 
written as 



u{a) 



1 



(3A- l)a 4 



tt 4 3r 2 -9s 

3 3 



(3.15) 



In Fig. [TJ we show the fate of the universe, which depends 
on the values of g Y and <? s . For the case of K = 1, there are 
two types of spacetime phases: One is as =>■ <BC, and the 
other is an oscillating universe. In fact, if g r > 0,g s < 
and g 2 + 12g s > 0, we find the scale factor a is bounded 
in a finite range as (0 <) <z m j n < a < a max (< 00), where 



g r - y/g* + 12g s 



g T + ^Jg 2 + 12g s 



(3.16) 



which gives an oscillating universe. The condition for an 
oscillating universe is written as 



9t > 



9£ 
12 



< g s < , 



(3.17) 
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which is defined by 



(a) K = 1 




(b) K 



FIG. 1: Phase diagram of spacetimes for A = 0. The os- 
cillating universe is found only for the case of K = 1. The 
stable and unstable static universes (s s and s u ) exist on the 
boundary Fo,i and Fo,-i, respectively. On Fo,-i, we also find 
dynamical universes with an asymptotically static spacetime; 

Su => »f, Su => 9f , B«=> Su, Or => Su- 



which is shown in Fig. \V[a) by "Oscillation" (the light- 
orange colored region) in the g s -g r plane. The equality in 
Eq. (|3.17|) . which is the curve r 0) i, gives a static universe 
with the scale factor a = as ■= \J <? r /6. 

For the case of K = —1, we find three types of space- 
time phases: <B<B => M (or M => bc), bb =>■ bc, and 
'Bounce (sec FigQJb)). On the boundary curve r .-i, 
which is defined by Eq. (|37T2|) . i.e., g s = 5l 2 /12 (g T < 0), 
we find an unstable static universe s u , and the dynamical 
universes with an asymptotically static spacetime; s u => 
bc, S u => M , ««=^> 5 U , or 5Vf s n - 

The bounce universe is found if g s < or g s = with 
g r < 0, which is shown by "Bounce" (the light-green re- 
gion). The radius at a turning point, ax is given by 



- {-g T + y/g? - 12ft 



(3.18) 



Next we shall evaluate the period of an oscillating uni- 
verse in the case of K = 1. The solution for Eq. p.4j) is 
given by 



da 



3A- 1 



E(<j)[a],k) , (3.19) 



E{4>,k):= [ d0Vl-k 2 sin 2 i 
fc and 0[a] are given by 



/ 2 2 
k . . V a max ~ a min 



a 2 - a 2 

"max " 



max mm 



The period T is given by 



T := 2(t n 



^max) " 2(2^ 



3A- 1 



(3.20) 

(3.21) 
(3.22) 

E(k), (3.23) 



where E(k) is the complete elliptic integral of the second 
kind defined by E(k) := E(ir/2, k). 

In order to evaluate the period, we consider some lim- 
iting cases, which are the boundaries of the region of 
Oscillation. In Fig. we show the potential u(a) by the 
blue curve for one boundary curve To,!, which is given 
by 5s = — 97/12. It gives a stable static universe with 
the scale factor as- We also show the potential near the 
other boundary of Oscillation (the positive (/ r -axis) by the 
dashed orange curve. Choosing, for example, g T = 1 and 
g s = —0.001, we find an oscillating univese with the scale 
factor a € [0.0316705,0.576481]. Since these two poten- 




FIG. 2: The potential «(5) for a stable static universe and 
an oscillating universe near the g r -axis. The "coupling" con- 
stants are chosen as g T = 1 and g s = —1/12 on To,! for a 
static universe, which radius is shown by as = 1/V&- We 
also show the case with g r — 1 and g s — —0.001 for an oscil- 
lating universe, which maximum and minimum radii are given 
by a max = 0.576481 and a m i n = 0.0316705, respectively. 



tials give the limiting cases, we find that < a m ; n < as 
and as < a max < v2as for an oscillating universe. 

In the limit of a static universe (near To,].), we find the 
period T5 as 




where E{<j>, k) is the elliptic integral of the second kind, 



(3.24) 

while in the other boundary limit (<7 S — > 0), we obtain 

(3.25) 
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Prom these evaluations, giving the value of g r , we find 
the period T of any oscillating universe is bounded in 
the range of (T ,T S ) for g s G (-.9 r 2 /12, 0). We then ap- 
proximate the period as T ~ gV 2 ■ 

We have found an oscillating FLRW universe because 
we have "negative" energy of "stiff matter" which comes 
from the higher curvature term. The condition for an 
oscillating universe is rewritten in terms of the original 
coupling constants as 



.93 + 3.92 > . 

?3 



(.93 + 3 52 ) 2 



< 9ffs + 3(7 6 + .97 < . 



(3.26) 
(3.27) 



B. A > (e = 1) 

In this case, the potential is given by 
1 



u(o) = 



(3A - l)a 4 



iv—4 ~6 9t~2 9s 

Ka —a a 



(3.28) 



For each value of K, we depict the fate of the universe in 
Fig |31 which depends on the values of <? r and g s . 



Su 



.9s 



0.4 

'Bounce 02 



*%, r l,l(+) 




Oscillation 
'Bounce 



Except for the case of <b<b =>■ ®c and a static universe, 
the expanding universe approaches de Sitter spacetime 
(exponentially expanding universe) because of a positive 
cosmological constant A. The oscillating universe exists 
if and only if K = 1 and the following conditions are 
satisfied: 



9r > 

gt 1]{ - ) (9r)<~9s 



< 

< 9s [1A]i + ) (~9r) 



(3.29) 
(3.30) 



where gj 1,1](±) is defined by Eq. with e = 1, K = 1. 

This condition gives the constraint on <? s as — 1/9 < g s < 
0. Note that in the limit of g Y <C 1 (i.e. A — >• 0), we 
recover the condition f|3.17|) . 

The boundaries of two different phases of spacctimcs 
consist of the <? r -axis, and two curves (rx^fi)) for K = 1 
or one curve (r"x _i(+)) for K = — 1. Those boundary 

curves T 1K ^ are defined by g s = (js (g T )- 

A stable static universe exist on the boundary curve 
while unstable static universes appear on the 
boundary curves r lj= |-x(-i-)- For K = 1, there are two 
types of static universes (stable and unstable) corre- 
sponding to two curves r 1 ! x(-) and rx,x(+), respectively, 
which coincide at g r = 1 and g s = —1/9. In the branches 
of unstable static universes (Ti^m)} we also find dy- 
namical universes with an asymptotically static space- 
time; S u => is or S u => «C on r lj x(+)(9s > 0); S n => is 
Or Su => bounce. s u on r ljl ( + )(.9 s < 0); s u => <£S or 5 U 
=*> «c on r lj _ 1 (+). 

The period T of an oscillating universe is calculated by 



T : = 2 



V-2«(5) 



(3.31) 
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(b)K = -l 

FIG. 3: Phase diagram of spacetimes for A > 0. The oscillat- 
ing universe is found only for the case of K = 1. The static 
universes (s s and s u ) exist on the boundaries T liK (±y We 
also find dynamical universes with an asymptotically static 
spacetime; s u => is or s u =>• ®C on r 11 ( + )(g s > 0); s u => is 
or Su => 'Bounce =>■ s u on r lil(+) (3 s < 0); Su => is or Su => %C 
on IY_i. 



where T = T/£, and flmax and flmin are the maximum 
and minimum radii of the oscillating universe. We shall 
evaluate the period near the boundaries of the parameter 
range of oscillating universes (the light-orange region in 
Fig. EJa)). We first show the potential u(a) for three 
(near-) boundary values of g s in Fig. U 

For the case with an unstable static universe (the 
dashed blue curve) (r lj x(+) with g s < 0), the larger dou- 
ble root of the equation of u(a) =0 is given by 



as 



:=^l(l + VT=t), (3.32) 



while the smaller root is 



CLT = CLrp 



- (l - 2 V ^ _ |7j , (3.33) 



We find non-singular evolution of the universe ('Bounce, 
Oscillation, or static) as well as the universe with a cosmo- 
logical singularity (a® =>■ <BC, as is, or ds =>■ ®c). 



which corresponds to a turning radius at a bounce. The 
period T diverges in the limit of a static universe, because 
flmax = 5s is the double root. 
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FIG. 4: The potential u(a) for a stable and unstable static 
universes (the solid blue and the dashed blue), and that for 
an oscillating universe near <? r -axis (orange). The constants 



are chosen as g T = 0.8 and g B = —0.0643206 on T 1 



-0.0245683 on T ± 



!(-)> 



and 



for static universes, which radii are 



given by as, and g T = 0.7 and g s — —0.001 for an oscillating 
universe, which maximum and minimum radii are given by 
5 max and a m i n , respectively. We also find s u =>■ (Bounce => s u , 
which bounce radius is given by a m i n . 



The period is then given by 



T 



3A- 1 



3A- 1 



x 2 sinh 1 

2 „■ 



l-(l-|ffr) S 



x < 



V3 



In 



9r 



2(l-| 5r ) 5 

(flr « 1) 



V3 



(3.39) 



(<7r 



The period To also changes from to oo along the g r -axis; 
~g s = (0 < g t < 3/4). 

1/2 

We summarize our result as T ~ g r when g r <C 1 , but 
it diverges near r lil( - + ), on which we have the unstable 
static universe. 



While, near a stable static universe (the solid blue 
curve) (r 11 (_)), the period is finite and is evaluated as 



T 8 = 



3A - 1 



3A- 1 



1/2 



X 7T 



1/2 



1/2 



l-(l-gr) 

3(1 -gr) 1 / 2 

r 5 y 2 



1/2 



(3.34) 



(9r « 1) 



V6 



The period T$ changes from to oo along the static curve 
r M(-)' 

The radius of this stable static universe is given by 
as = , which is the smaller root of the equation 

of u{a) = 0. The larger root clt = S,g' 1 ^ + ' ) corresponds 
to a turning radius of a bounce universe, which is shown 
by or in Fig. [4] 

There is another boundary limit, i.e., g s — > CP. In this 
limit, we find the roots of <u (a) = as 



~2 







r 1 

















In 





(3.35) 
(3.36) 

(3.37) 



Since the largest root (0,3) corresponds to a turning radius 
cit of a bounce universe, the oscillation range is [01,02], 
and then the period is evaluated approximately by 



Tn = 2 



do 



y/-2u{a) 



(3.38) 



C. A<0(e = -1) 
The potential is given by 
a (a) 



1 



(3A - l)a 4 
We summarize our result in Fig. [5] 



Ka + a a 

3 3 



, (3.40) 









m - 


.0 


05 




-0B2 




No Sol 




0.06 



Oscillation. 




(a) K = 1 



Oscillation 
<B<B ^ ®C 



<^./, 


gs 


, <m 'Be 


i(+) 




-1.0 " -05 




05 1 




. -02 


Oscillation 


No Solution 








-0.6 


\ r -l,-l(-) 




\5s 



(b)K = -l 

FIG. 5: Phase diagram of spacetimes for A < 0. The os- 
cillating universe is found for both K = ±1. The static 
universe exists on the boundary r_i j i(_) (K = 1) and on 
r_ 1 ,_ 1 (±) [K — —1). In the branch of unstable static uni- 
verse on r_i if_|_), we also find dynamical universes with an 
asymptotically static spacetime; SB =>• s u , S u => "BC, or 5 U => 

Bounced Su- 
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In this case, if g s > 0, we find a big bag and a big 
crunch singularities («S => $c) except for a small re- 
gion in K = — 1. If g s < 0, however, we always find an 
oscillating universe if the solution exists. 

The conditions for an oscillating universe is shown by 
the light-orange region in Fig. [5j which is given by the 
following inequalities: 
For K = 1, 

g r > 

fli- 1 ' 1]( - ) (ffr)<5.<0, (3.41) 

and for K = —1, 

g[-b-l](-)(g r ) < g s < with g r > 0; 

9!" lrl|H (ft) < 5s < with g r < 0. 

(3.42) 

In the limit of g r <C 1 (i.e. A —> 0) for K = 1, we recover 
the condition (|3.17l) . 

The boundary of the range of oscillating universe is 
given by the positive g r -axis, and r_x .i(-) for K = 1, and 
r_i ,-i(±) for K = —1. On those boundaries r_x .k(±), 
which are defined by g s ~ g[ 1 \<?r) [K = 1) and g s = 

gl 1: (<?r) (K = —1), we find a stable and unstable 
static universes. 

The period of an oscillating universe is given by 
Eq. p.31[) . We again evaluate its value near the boundary 
curves (T_i ,k(-)) an d the positive <7 r -axis. The poten- 
tials u (a) for the (near-) boundary values of g s are shown 
in Fig. E] (if = 1), and Figs. and [8] (K = -1). 




The turning point is given by or = 62, where 

a 2 2 + 2^/T+g\) . (3.44) 




FIG. 7: The potential 11(a) for a stable static universe (blue) 
and an oscillating universe near y r -axis (orange) for K = — 1. 
We set g-c = 0.2 and g B = —0.581008 on r_i,_i for a static 
universe, which radius is given by ds = 0.835752, and g T = 0.2 
and g s — —0.001 for an oscillating universe, which maximum 
and minimum radii are given by a max = 1.03075 and a m i n = 
0.0683656, respectively. 




-0.4 



FIG. 8: The potential 11(a) for a stable and unstable static 
universes (blue and red, respectively) and an oscillating uni- 
verse on <? r -axis (dashed orange) for K = — 1. We set g r = 
-0.5 and g s = -0.134123 on r_ 1 _ 1( _ ) , and g s = 0.0230119 
on r_i i _i( + ) for static universes, which radius is given by 
as = 0.754344, and g r = —0.5 and g s — for an oscillat- 
ing universe, which maximum and minimum radii are given 
by a max = 0.888074 and a min = 0.459701, respectively. We 
also find Su =>• 'Bounce => Su, which bounce radius is given by 
a T = 0.897072. 



FIG. 6: The potential 11(a) for a stable static universe (blue) 
and an oscillating universe near <? r -axis (orange) in the case 
of K = 1. We set g r = 0.8 and g s = -0.0477674 on r_ M for 
a static universe with the radius 5s = 0.337461, and g T — 0.8 
and g s = —0.001 for an oscillating universe, which maximum 
and minimum radii are given by 2 max = 0.466615 and 5 m i n = 
0.035439, respectively. 



Note that the period diverges in the limit of an un- 
stable static universe (on r_i .-i(+)), where we find the 
radius of a static universe by as = Si 



1 



Near a stable static universe (F_x i i(_) and r_x _i(_)), 
the period is evaluated as 



,'3A- 1\ 1/2 



3(l+5r)V2 

which approaches a constant 

/o\ 1 \ 1/2 
7T / 3A — 1 N 



V3 V 2 



(3.43) when g T > 1. 



1/2 



(3.45) 



(3.46) 
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Near the lower bound of g r , we find 



gl/2 

7T /3A - 1 \ 1/2 I -> (as g T -> for K = 1) 



(1 + .9r) _ 7 -> oo (as g r -> -1 for A" = -1) . 



Hence the period Tg changes from to a finite value 
(|3.46p along the curve r_ Xi i(_) for K = 1 , while from oo 
to at the same finite value along the curve r_ li _ 1 /). 
The radius of a static universe is given by 



\(VT+h-K). (3. 



- _~[-l,K]{-) 
a s — a s 



In the case of g r < — 3/4 with A = — 1 , there is another 
zero point of u (a), which gives a maximum turning point 
of ®3 «c, i.e., 



or 



gtr 1 .- 1 "-) = ^ (1-271H). (3.49) 




Near g r -axis, we find the solutions of the equation 
u (a) = as 



(3.50) 



as well as do ~ 0. Wc have a maximum radius 5 max = d + , 
and find that the minimum radius a m j n is almost zero for 
g T > because a 2 _ < 0, but in the case of A' = — f, 
for —3/4 < 5 r < 0, we find a finite minimum radius 

"min — GE— ■ 

Using those values, we evaluate the period as 



/ 3A- 1 



sec 



4 

1 + 3* 



(3.51) 



for K = 1, and 
/3A- 1 



7T — sec 



1 + jjft (ft > 0) 

(-3/4 < ~g T < 0) 



for A = — 1. The period To also changes from to 00 
along the g r -axis; g s = (0 < g T < 3/4). 

In the case with the detailed balance condition, since 
A < 0,g T = -9/4, g s = 0, wc do not find any FLRW so- 
lution. If wc include matter fluid, the result will change. 
For example, if we have "radiation" fluid, which energy 
density is proportional to a~ A , we should shift the value 
of _o r . Then if —3/4 < p r < 0, we find an oscillating uni- 
verse for A = -1, which period is tt[(3A - l)/2] 1 / 2 . The 
equality (g r = —3/4) gives a static universe. 



IV. TOWARD MORE REALISTIC 
COSMOLOGICAL MODEL 

In the Hofava-Lifshitz gravity without the detailed bal- 
ance condition, we find a variety of phase structures of 
vacuum spacetimes depending on the coupling constants 
g r and g s as well as the spatial curvature A and a cosmo- 
logical constant A. Note that there is no vacuum FLRW 
solution in the case with the detailed balance condition. 
We summarize our result in Table |U We have obtained 
an oscillating spacetime as well as a bounce universe for a 
wide range of coupling constants. We have also evaluated 
the period of the oscillating universe. 





K = 1 


K 


-1 


A > 


* OsciUation 

* dS'i VRouncc 

* <B1 =^-BC 

* 23 => dS (dS => <BC) 


* £S<*=>S ounce 

* mm => ds (ds => mc) 


r l,l(±) 


* Su, S s 
I (Su => 9C) 

iS (<«=>• Su) 
> 'Bounce 


r i--i(+) 


* Su 

(s u =>■ mc) 


* mm^s 

* Su => 

* Su ^= 


* mm=>Su 

* Su =>• ds 


A = 


* OsciUation 

* B<8 =>SC 


* M <= 

* mm =>sc 

* mm => M (m => mc) 


r .i 




r o,-i 


1 (in => Be) 

(stf 5u) 


* Ss 


* mm => Si 

* Su =>• M 


A < 


# OsciUation 

* <B<B =5-mC 


* OsciUation 

* mm =>2sc 


r -l,l(-) 




r -i,-i(±> 


i (s u => mc) 

(Bounce 


* Ss 


* mm => 5i 

* Su •<=>■ 



TABLE I: Summary: What type of spacetime is possible for 
each A and each K. Non-singular universes are shown by the 
colored letters (an oscillating universe and dynamical space- 
times evolving in a finite scale range by red, static universes 
by blue, dynamical spacetimes evolving from or to an asymp- 
totically infinite scale by green), ds, %'B, <BC, Su, S s and M 
denote de Sitter space, a big bang, a big crunch, an unstable 
static universe, a stable static universe, and Milne universe, 
respectively. 



In our analysis, we assume that the integration constant 
C from the projectability condition vanishes. If C ^ 0, 
one may find a different story. In fact, if <7 S = and g r < 
just as the case with the detailed balance condition, we 
will find the similar vacuum solutions to the present ones, 
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because C and g Y without g s -term play the similar roles 
to those of g T and g s in the present model. For example, 
we obtain an oscillating universe for large C(> 0) with 
g s = 0, g r < 0, A = and K = 1. This avoidance of 
a singularity is, however, caused by the negative "radi- 
ation" density from the higher curvature terms. Hence 
if one includes the conventional radiation, then the ef- 
fective g T becomes positive as we will show below, and 
as a result the universe will inevitably collapse to a big- 
crunch singularity. Furthermore, if radiation field evolves 
as a~ 6 in the UV limit |23|. the inclusion of such radiation 
will kill the possibility of singularity avoidance by "dark" 
radiation. 

As we have evaluated, the oscillation period and am- 
plitude are expected to be the Planck scale or the scale 
£ defined by a cosmological constant A, unless the cou- 
pling constants arc unnaturally large. Hence it cannot 
be a cyclic universe, which period is macroscopic such as 
the age of the universe. 

In order to find more realistic universe, we have to 
include some other components, which we shall discuss 
here. First of all, one may claim inclusion of matter fluid. 
When we include a dust fluid (P = 0), the conventional 
radiation (P = p/3), and stiff matter (P = p), we can 
treat such a case just by replacing the constant g^, g r 
and g s with 



ZL 



fid 
<h 



8C + <?dust 

6(ff3 + 3.92) +5rad 

12(9 ff5 + 3.9 6 + g 7 )K 



ffstiff 



(4.1) 



where gdust, ffrad and .9 s tiff, which come from real dust 
fluid, radiation and stiff matter, are positive constants. 
In this case, the present analysis is still valid. If .9 ra d 
is large enough just as our universe, a maximum scalar 
factor a max of the the oscillating universe will become 
large (see, for example, Eq. (|3.16|l ). and then it can be a 
cyclic universe. 

If the equation of state is still given by P = wp 
(inconstant), the analysis is straightforward. When we 
have other types of matter fields, e.g. a scalar field with 
a potential, the analysis will be more complicated. The 
phase space analysis may be appropriate for the case with 
a scalar field [63 1. 

From our present analysis, one may speculate the fol- 
lowing "realistic" scenario for the early stage of the uni- 
verse. Suppose a closed universe is created from "noth- 
ing" initially in an oscillating phase (see Fig. [9]) [HI, [6|| . 
Such a universe may be very small and oscillating be- 
tween two radii (a m ; n and a max ) with a time scale £. If 
we have a positive cosmological constant (A > 0), there 
exists a potential barrier as shown in Fig. [9] 

After numbers of oscillations, the universe may quan- 
tum mechanically tunnel to a bounce point ax- Then 
the universe will expand to de Sitter phase because a 
positive cosmological constant, finding the universe in a 

















a 








1.0 




/ a™, 







FIG. 9: 



macroscopic scale 1 . Furthermore, one can refine this sce- 
nario, if there exists a scalar field, which is responsible 
for inflation, instead of a cosmological constant. Before 
tunneling, we may find the similar scenario to the above 
one. After tunneling, the potential of the scalar field 
will behaves as a cosmological constant in a slow-rolling 
period. We will find an exponential expansion of the uni- 
verse after tunneling. However, inflation will eventually 
end and the energy of the scalar field is converted to that 
of conventional matter fluid via a reheating of the uni- 
verse. We find a big bang universe. Since the universe 
is closed, but the scale factor has lower bound because 
of negative "stiff matter" , we will find a macroscopically 
large cyclic universe after all. To confirm such a scenario, 
we should analyze the dynamics of the universe with an 
inflaton field in detail. The work is in progress. 

We also have another extension of the present FLRW 
spacetime to anisotropic one. It may be interesting and 
important not only to study the dynamics of Bianchi 
spacetime (6(| [67j but also to analyze the stability of the 
FLRW universe against anisotropic perturbations [68]. 
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Appendix A: stability of a flat background and the 
coupling constants 

In this Appendix, we discuss the conditions on the cou- 
pling constants by which gravitons are perturbatively sta- 



1 After we have written up this paper, we have found [60H , in which 
a cosmological transition scenario from a static (or an oscillating) 
universe to an inflationary stage was discussed. They assume 
that the equation of state changes in time, which mechanism is 
not specified. 
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blc. From the perturbation analysis around a flat back- 
ground, we obtain the dispersion relation for the usual 
helicity-2 polarizations of the graviton (l7j . 



J TT(±) 



= -9i k 2 



93 



M 2 



±34 



99 



AT 4 



(Al) 



The stability both in the IR and UV regimes requires 



gi < 0, .99 > 0. 



(A2) 



By a suitable rescaling of time, we then set g± = — 1. 

As a result of the reduced symmetry (|2.4j) the longitu- 
dinal degree of freedom of the graviton appears, and its 
stability is more subtle. First of all the longitudinal gravi- 
ton is plagued with ghost instabilities forl/3<A<l[l[. 
The dispersion relation for the longitudinal mode turns 
out to be [13 



3A- 1 
A - 1 



w£ = gik 2 + (832 + 393) 



M 2 



+(- 



339 



k 



M 4 



(A3) 



We see that the sound speed squared is negative in the IR 
if gi < and A > 1, which implies that the longitudinal 
graviton is unstable in the IR |36[. However, this fact 



itself does not necessarily mean that the theory suffers 
from pathologies, because whether or not an instab ility 
really causes a trouble depends upon its time scale [271 ]. 
Moreover, there is an attempt to improve the behavior 
of the longitudinal graviton by promoting N to an in- 
dependent function and adding terms constructed from 
the 3-vector diN/N in the Lagrangian (35|. 2 It can be 
shown that the non-pro jectable Hofava gravity thus ex- 
tended appropriately does not plagued with instabilities 
of the longitudinal gravitons [35||. In light of these sub- 
tleties, we do not consider the stability of the longitudinal 
sector furthermore, while we do require the stability for 
the usual helicity-2 polarizations of the graviton. 

Note that the detailed balance condition satisfies g\ < 
and gg > 0. 



Appendix B: quantum tunneling from an oscillating 
universe 

In the case of K = 1 and A > 0, we have a bouncing 
universe as well as an oscillating universe. These two so- 
lutions are separated by a finite potential wall as we see 
in Fig [5J Hence we expect quantum tunneling from an 
oscillating universe to an exponentially expanding uni- 
verse. In this Appendix, we shall evaluate the tunneling 
probability. 



2 Obviously, in this case the Hamiltonian constraint is imposed 
locally and the additional dust-like component does not appear 
in the Friedmann equation. 



First we consider the normalized Euclidean metric 



dS z = df z + Y{f)dH z K=1 
which satisfies the following equation 
V 2 -2u{b) = 0, 



(Bl) 



(B2) 



where the prime denotes the derivative with respect to 
the Euclidean time f , and the potential 11 is written as 



2u(b) = 



2 1 
3A- 25J 



-(b 2 -b 2 m ^)(b 2 -bl in )(b 2 -b 2 T ) (B3) 



The variables with a tilde arc normalized ones by use of 
the scale length I = \/3/A just as in the text. The 
bounce solution 6(f) is obtained by integraton of Eq. 
(|B2|) . The Euclidean action is given by 



(B4) 



S E = 3(3A - 1)£ J dfd 3 xb ^b' 2 + 0(6) 
Using Eq. (|B2[) , we find the action Se as 



S E = 3(3A - l)tV 3 / dbbJ2u(b) 



(B5) 



where V3 = 2ir 2 is the volume of a unit three sphere. 
Introducing u by 



b 2 = b 2 T (l - k 2 u 2 ) , 
where k 2 = (b 2 ^ - 6 2 n ax)/^T(< !)• We then find 

127r 2^2 _ 

S E = -0 ( b T ~ frmax) 2 ^ ~ fo min) 



(B6) 



1 u 2 du 
1 — k 2 u 2 



- m 2 u 2 ) , (B7) 



where m 2 = (b\ - b 2 a&x )/(b 2 T - &^J(< 1). 

It can be easily evaluated in the limit of a static uni- 
verse, i.e., g s = 3l 1,11(_:) (.9r)- Using 6 max sa 6 min w b s , we 
find 

b E - 5 — (b T - b s ) 



3 -2k 2 1 - k 2 



3fc 4 



k 5 



tanh k 



(B8) 



where k = Jb% - b 2 /b T . Since &| = (1 + 2^1 - ff r )/3 
and b\ — b 2 s = \J\ — g T , we find 



S E = ^(Uj r ) 1/4 



(1 + 2V1— Yr) 1/2 (i-VT^ j,) 

V3(l-g r )!/4 



tanh 1 k 



(B9) 
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with 



■VI - ffr 

1 + 2VT - fx 



(BIO) 



If the vacuum energy (or potential) just after tunnel- 
ing is the Planck scale, the probability is evaluated as 



P 



-(60-120) 



, which is very small but finite. 



The tunneling probability is given by P ~ e Se . 
We show the behavior of Se in Fig. [KjJ We find 



P 



exp 



exp 



-(20-40)x j- 



(60 - 120) x 



'PL 



(Bll) 



except for two limiting cases: g T ~ 1, in which Se van- 
ishes, and g r ~ 0, in which Se diverges. In the for- 
mer case, the potential barrier vanishes giving a high 
tunneling probability, while in the latter case, the po- 
tential barrier diverges giving zero tunneling probability. 



i-3 



Co 20 



9r 



FIG. 10: 
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